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ABSTRACT 

We present a finely-binned tomographic weak lensing analysis of the Canada-France- 
Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from 
the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model 
and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, 
with accurate shear and photometric redshifts for a galaxy sample with a median redshift of 
z m = 0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six 
redshift bins, each spanning the angular range of 1.5 < 9 < 35 arcmin. We combine this 
CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with 
data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the 
Hubble constant from the HST distance ladder. This leads to constraints on the normalisation 
of the matter power spectrum er 8 = 0.799 ± 0.015 and the matter density parameter Sl m = 
0.271 ± 0.010 for a flat ACDM cosmology. For a flat wCDM cosmology we constrain the 
dark energy equation of state parameter w = —1.02 ± 0.09. We also provide constraints for 
curved ACDM and uiCDM cosmologies. We find the intrinsic alignment contamination to 
be galaxy-type dependent with a significant intrinsic alignment signal found for early-type 
galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is 
found to be consistent with zero. 
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1 INTRODUCTION 

Cosmological weak lensing is the study of the weak gravitational 
distortions imprinted on the images of distant galaxies by large- 
scale structures. A series of deflections are induced by the gravita- 
tional potential of the dark and luminous matter that light passes, as 
it travels through the Universe. This lensing effect results in a co- 
herent distortion, detected in the observed images of galaxies, that 
allows us to infer the distribution and density of matter in the Uni- 
verse. This well understood physical effect is recognized as one of 
the most powerful probes of cosmology, allowing not only the di- 
rect study of dark matter, but also, through the study of the growth 
of structures, a unique probe of gravity and dark energy on large 
scales (see review by Weinberg et al. 12012] and references therein 
for more details). 

It has long been recognized that the optimal way to extract 
cosmological information from the detection of weak gravitational 
lensing is to utilize redshift information, for example by separat- 
ing the lensed galaxy sample into a number of tomographic bins 
using photometric redshift information (Hu 19991. This idea was 
explored theoretically to determine the optimal number of redshift 
bins (see for example Huterer 2002; Simon et al. 2004) and also 
significantly extended to consider a fully three-dimensional analy- 
sis ((Heavens 2003). These early predictions suggested that it was 
possible to achieve up to an order-of-magnitude improvement on 
cosmological constraints when weak lensing measurements were 
made in combination with photometric redshifts. This drove weak 
lensing survey designs to include multi-band optical imaging for 
redshift estimates, in addition to high-resolution imaging in a sin- 
gle band for the measurement of the weak lensing distortions. 

In order to detect cosmological weak lensing, correlations are 
measured between the shapes of galaxies whose observed elliptic - 
ity, e bs, is related to their intrinsic ellipticity, e s , and the weak cos- 
mological shear distortion that we wish to extract, 7, througrQ 

£obs = e s + 7 • (1) 

The observed angular two-point correlation function (eobs e obs) is 
then determined by averaging over all galaxy pairs (a, b) separated 
by angle 9. This observed quantity is related to the angular two- 
point shear correlation function (7®7 ) through 

(eSbs4bs) = <es£s) + (£sj b } + (l a 4) + <7 a 7 6 > ■ ( 2 ) 

For the majority of weak lensing analyses to date, the first three 
terms on the right hand side of equation [2] have been assumed to 
be sufficiently small to be ignored (although see discussion below 
for analyses that do not make this assumption). In this case the ob- 
served angular two-point correlation function (£oba e obs) i s equated 
with the cosmological shear correlation function (7 a 7 6 ) , which can 
be directly related to the underlying matter power spectrum of den- 
sity fluctuations in the Universe (see Section [3] for more details). 
It is, however, possible that these intrinsic terms are significantly 
non-zero, arising from correlations induced during galaxy forma- 
tion, between a galaxy's intrinsic shape and its local density field. 
It was noted early on that such an effect would be significantly 
more detrimental to the future measurement of weak lensing in to- 



1 For clarity we assume this simplified form to relate galaxy ellipticity and 
shear. In detail, the relationship also depends on the ellipticity estimator 
used in the analysis and how rapidly the galaxy ellipticity varies at different 
isophotal limits (see Section |3~T} . We also assume here the weak shear limit 
that |7| <C 1. 



mographic bins ( King & Schneider 2003 ; Heymans et al. 2004 1, in 
comparison to the standard two-dimensional analysis of that time. 

The study of the impact of intrinsic galaxy alignments on 
weak lensing studies initially focussed only on the intrinsic align- 
ment of physically nearby galaxies; (£563), hereafter referred to as 
TP . The broad agreement in results between the first estimates from 
numerical simulations ( Heavens et al. 2000 , Croft & Metzler 2000 1, 
analytical studies (Catelan et al. 2001 ; Crittenden et al. 2001 ; Lee 
|& Pen||2001| >, and the first low redshift observational constraints 
( |Pen et aLp OOO , Brown et al. 2002 ) resulted in a consistent picture; 
for deep weak lensing surveys, the contamination to the weak lens- 
ing signal was expected to be less than a few per cent effect. [Hirata| 
& Seljak (2004) were the first to highlight, however, the impor- 
tance of also including the shear-shape correlations in the analysis, 
(e s 7), for galaxies that are separated by large physical distances 
along the line of sight. In this case the background galaxy experi- 
ences a shear 7 caused by the foreground tidal gravitational field. 
If the foreground galaxy has an intrinsic ellipticity that is linearly 
correlated with this field, (e s 7) is no longer expected to be zero. 
We refer to this effect hereafter as 'GP. 

The most recent observational results have focused on infer- 
ring the amplitude of the II and GI signal by measuring the local 
cross-correlation between galaxy number densities and ellipticities, 
to determine the correlation between galaxy shape and the local 
density field. This method was first implemented using low red- 
shift z ~ 0.1 galaxies in the Sloan Digital Sky Survey spectro- 
scopic sample (SDSS, Mandelb aum et af1 [2006). It was then ex- 
tended to higher redshifts using the SDSS Luminous Red Galaxy 
sample (LRG) between 0.15 < z < 0.3 5 jHirata et al.|2007), and 
the MegaZ-LRG sample out to z ~ 0.6 (loachi mi et al.|2011| >. In 
these extensive studies, a clear dependence on galaxy type is de- 
tected, with the most massive red galaxies exhibiting the strongest 
intrinsic correlations. Hirata et al. (2007 ) show their data match 
models for this galaxy-type dependence as predicted from numer- 
ical simulations (Hey mans et al.|20 06). The massive LRG popula- 
tion is, however, in the minority when it comes to a typical deep 
cosmological weak lensing survey, which is dominated by the blue 
galaxy population. The most representative observational analysis 
of this effect therefore comes from the highest redshift measure- 
ments of the GI and II effect, using blue galaxies from the Wig- 
gleZ survey out to z ~ 0.7 (Mandelbaum et al. 2011). The re- 
ported null detection of this effect for blue galaxies results in a 
predicted systematic error of at most ±0.03 on the amplitude of 
the matter power spectrum, <rg, for a non- tomographic analysis of 
a survey like the Canada-France-Hawaii Telescope Lensing Survey 
(CFHTLenS). We perform a tomographic weak lensing analysis of 
this survey in this paper. 

For a deep survey, such as CFHTLenS, the correlation be- 
tween a galaxy's shape and its local density field impacts more 
strongly on a tomographic analysis than a two-dimensional anal- 
ysis for the following reasons. For a deep two-dimensional analy- 
sis, there is only a small fraction of galaxy pairs, at fixed angular 
separation 8, that are physically close enough to have formed in 
the same density field and hence experience some degree of align- 
ment. The significant II signal from these close pairs is therefore 
greatly diluted as the majority of galaxy pairs are well separated 
in three-dimensions. The lensing of background galaxies by a fore- 
ground structure is most efficient when the foreground structure 
is well separated from the source. In a two-dimensional analysis, 
the signal from the galaxy pairs that experience the strongest GI 
effects is therefore again diluted by the presence of closer galaxy 
pairs in the analysis where the lensing effect is less efficient. In 
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contrast, a tomographic analysis will include redshift bin combina- 
tions that, in fact, enhance the contamination by both these effects. 
The auto-correlation measurement of the lensing signal within nar- 
row redshift bins dramatically increases the fraction of physically 
close pairs and hence the II contribution. The cross-correlation of 
well separated redshift bins, where the lensing is most efficient, en- 
hances the GI contribution. For this reason we cannot ignore the 
intrinsic ellipticity terms in equation [2] when undertaking a tomo- 
graphic weak lensing analysis and need to apply methods to mit- 
igate the impact of this astrophysical contamination to the signal. 
These mitigation strategies unfortunately limit the dramatic order- 
of-magnitude improvements initially anticipated from tomographic 
weak lensing analyses, but we still expect to find improved con- 
straining power when implementing a tomographic analysis {Bridle| 
|& King|2007) . 

A series of different proposals to mitigate the impact of intrin- 
sic alignments on weak lensing measurements have been made in 
the literature, and in some cases applied to data. An optimal weight- 
ing scheme to down-weight physically close pairs was proposed by 
Heymans & Heavens ( 2003 and applied to a weak lensing analysis 



of the COMBO- 17 survey ( |Heymans et aTp 004 see also King & 



Schneider 2002). This method, however, only negates the II con- 
tamination and requires prior knowledge of the angular and red- 
shift dependence of the intrinsic ellipticity correlation function in 
order to optimally analyse the data. Motivated by our lack of knowl- 
edge in this area, Joachimi & Schneider ( 2008 ) proposed a nulling 
method where only the characteristic redshift dependence of the II 
and GI correlations in different combinations of redshift bins (see 
discussion above) are used to derive an optimal weighting scheme 
that, in an ideal case, nulls all II and GI contributions to the to- 
mographic analysis. This method, however, has been shown to sig- 
nificantly degrade cosmological parameter constraints (Joachimi & 
|Schneider|2009l >. 

The main alternative to mitigating intrinsic alignment contam- 
ination by using different weighting schemes, is instead to use a 
model fitting approach, and it is this technique that we exploit in 
this paper. This approach was first highlighted by King ( 2005 1 who 
demonstrated that a simultaneous model fitting analysis of finely 
binned tomographic data, using a sufficiently flexible parametrized 
model for the II and GI signals, allows for the marginalisation over 
the II and GI nuisance parameters in the final cosmological analy- 
sis. This idea was extended by Bridle & King ( 2007 ) to determine 
how the figure of merit for future cosmological surveys would de- 
grade based on the flexibility of the II and GI model, or how much 
prior knowledge we are prepared to assume, and the accuracy of 
the photometric redshifts for the survey. For a typical photometric 
redshift error of a z ~ 0.05(1 + 2;) the figure of merit for measuring 
dark energy would decrease by 20 to 50 per cent depending on the 
allowed flexibility of the model. Similar conclusions were drawn 
by |Kitching et al.| p008} and |Kitching & Taylor| pOTT) when in- 
vestigating the loss of constraining power when intrinsic alignment 
modeling was marginalised over in a tomographic or fully three- 
dimensional weak lensing analysis, and more recently by Blazek 
|et al.| f2012) who investigated the impact of marginalising over in- 
trinsic alignments in a galaxy-galaxy lensing analysis. The model 
fitting approach does, however, become more promising if the lens- 
ing data are analysed simultaneously with galaxy clustering data. 
These extra data act to self-calibrate the II and GI signals (Zhang 
[20T0{ | Joachimi & Bridle|2^1^ . [KlrFeTar^ ( [2^T0) present the first 
example of such a joint analysis, combining two-dimensional weak 
lensing data from the 100 square degree lensing survey ( Benjamin 
|et al.|[2007} and SDSS shear-shape clustering data (Mandelbaum 



|et al. 12006} . A single parameter analytical model for intrinsic align- 
ments from |Hirata & S eljak (2004) is used in this analysis, and we 
describe this 'non-linear intrinsic alignment' model in more detail 
in Section |X2| Huff et al. (201 1) also use this model to remove GI 
contamination to their cosmic shear analysis of the SDSS-Stripe 
82 survey, but they fix the contamination using a mean measure- 
ment of the total amplitude for the GI signal as determined by |Hi-| 
|rata et al.| l |2007^ . Finally, [Fuet al. 1(2008) also mitigate the impact 
of GI on their cosmic shear analysis of the third year data from 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS), pre- 
senting a simultaneous two-dimensional weak lensing and intrinsic 
alignments model fitting method. This analysis however uses an al- 
ternative model for the GI contamination, motivated by numerical 
simulations (Heymans et al.|2006| >. 

A third way to account for the intrinsic alignment signals 
has been proposed, modifying the covariance matrix such that a 
marginalisation over possible functional forms of the II and GI 
power spectrum is permitted (Kitching & Taylor 2011). This has 
the advantage in that all functional forms are explored, and does not 
require explicit estimation of any nuisance parameters, but makes 
the assumption that the variance of the intrinsic alignment functions 
is Gaussian. 

The first combined weak lensing and photometric redshift 
analyses to directly detect the growth of structure, came from the 
COMB O-l 7 survey ( |Bacon et al.|2005||Krtching et al.|2007) and 
CFHTLS-Deep fields ( Semboloni et al. 2006 ). These analyses were 
followed by two independent tomographic analyses of the Hubble 
Space Telescope COSMOS survey (Massey et al.|2007||Schrabback| 
|et al.|2010| i. The areas of all these surveys were considered suffi- 
ciently small, all less than a few square degrees, such that any con- 
tributions from intrinsic alignments could be ignored in comparison 
to the large statistical errors. |Schrabback eTal] {2010} did, how- 
ever, attempt to mitigate any errors by removing all auto-correlated 
narrow redshift bins from their analysis to reduce the impact of 
the enhanced II signal in those bins. In addition they removed all 
luminous red galaxies from their galaxy sample as the shapes of 
this type of galaxy have been shown to be the most strongly corre- 
lated with the local density field ( Joa chimi et al.|2011) . As survey 
sizes grow and statistical errors decrease, it is not possible to ig- 
nore intrinsic alignments when analysing weak lensing in redshift 
bins. The recent two-bin tomographic analysis of the 154 square 
degree CFHTLenS, prese nted in|Benjamin et al. (|2013|>, therefore 
combines the strategies of |Huff et al.| ( |201 l) and Schrabback et al. 
( |2010} to mitigate intrinsic alignment contamination. Using the lin- 
ear tidal field intrinsic alignment model of |Hirata & Se ljak (20041, 
and following Bridle & King ( 2007 ) by fixing its amplitude to the 
observational constraints obtained by |Brown et al. ( 2002), they esti- 
mate the II and GI contamination to the cosmic shear measurement. 
They then limit their analysis to two broad high redshift bins with 
photometric redshifts 0.5 < z p h < 0.85 and 0.85 < z p h < 1.3 
such that any contamination from intrinsic alignments is expected 
to be no more than a few per cent. This broad bin tomography mea- 
surement is also used to constrain parametrized modified gravity 
models in Simp son et al.| (f20 1 3 >. In this paper we use the same 
CFHTLenS data, presenting the first tomographic weak lensing 
analysis to apply a full model fitting approach to mitigate the im- 
pact of intrinsic alignment contamination on shear correlation func- 
tions. A fully three-dimensional weak lensing analysis (Kitching 
is applied to the same CFHTLenS data in Kitching 



et al.| 


2011 


et al.|( 


2013 



This paper is set out as follows. In Section [2] we describe 
CFHTLenS and the auxiliary data sets used in this analysis. We 
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outline our methodology and chosen intrinsic alignment model in 
Section|3] additionally describing how our tomographic analysis is 
constrained by our requirements on the accuracy of the covariance 
matrix estimated from N-body lensing simulations. We present our 
results in Section [4] comparing joint parameter constraints from 
different combinations of CFHTLenS data with the cosmic mi- 
crowave background data from WMAP7, baryon acoustic oscilla- 
tions data from BOSS, and a prior on the Hubble constant from the 
HST distance ladder. In Section|5]we focus on the constraints that 
can be placed on the amplitude of the intrinsic alignment signal 
for early-type and late-type galaxies with this type of cosmological 
parameter analysis, with concluding remarks in Section|6] 



2 THE CANADA-FRANCE-HAWAII TELESCOPE 
LENSING SURVEY 

The Canada-France-Hawaii Telescope Lensing Survey 
(CFHTLenS) is a 154 square degree deep multi-colour u* g'r'i'z' 
survey optimised for weak lensing analyses, observed as part 
of the CFHT Legacy Survey (CFHTLS) on the 3.6m Canada- 
France-Hawaii telescope. The data span four distinct contiguous 
fields: Wl (~ 63.8 square degrees), W2 (~ 22.6 square degrees), 
W3 (~ 44.2 square degrees) and W4 (~ 23.3 square degrees). 
The CFHTLenS analysis of these data presents the current state- 
of-the-art in weak lensing data processing with THELI (Erben 
|et al.||2013} , shear measurement with lens/rf ( |Miller et al.||2013} 7 
photometric redshift measurement from PSF-matched photometry 
(Hildebrandt et al. 20121 using the Bayesian photometric redshift 
code BPZ iBenftez 2000), and a stringent systematic error analysis 
( |Heymans et al.|2012[ >. The resulting galaxy catalogue that we use 
in this analysis includes a shear measurement e b s with an inverse 
variance weight w and a photometric redshift estimate zbpz with a 
probability distribution P(z) and best-fit photometric galaxy type 
Tbpz- We apply the galaxy size and signal-to-noise dependent 
shear calibration corrections described in Miller et al. ( 2013} and 
Heym ans et al.| ([2012), and only use the subset of 75 per cent of 
the survey data that have been verified as science-ready and free of 
significant systematic errors. This has been demonstrated through 
a series of rigorous cosmology-insensitive tests on both the shear 
and photometric redshifts measurements, in combination (see 
|Heymans et al.|2012| for the full details). |Benjamin et al.| {2013} 
also use a cross-correlation analysis to verify the accuracy of the 
measured redshift distributions P(z) when the galaxy sample is 
limited to those galaxies with a most probable photometric redshift 
estimate between 0.2 < zbpz < 1-3. In light of these analyses, 
that demonstrate the robustness of these data to systematic errors, 
we do not present any further systematic error analyses in this 
work, referring the reader to Heymans et al. (2012), and references 
therein. For the redshift selection 0.2 < zbpz < 1.3, the galaxy 
sample has a weighted mean redshift of z — 0.75, and a weighted 
median redshift of z m = 0.70, as determined from the weighted 
sum of the P(z). The effective weighted galaxy number density, in 
this redshift range, is n c s = 11 galaxies per square arcmin. 



2.1 Auxiliary cosmological data 

In this analysis we present joint cosmological parameter con- 
straints by combining our tomographic weak lensing analysis of 
CFHTLenS with up to three complementary data sets to break pa- 
rameter degeneracies. We include the temperature and temperature- 
polarisation cosmic microwave background power spectra from the 



Wilkinson Microwave Anisotropy Probe ^Larson et al.|201 l) here- 
after referred to as WMAP7). We incorporate the measurement of 
baryonic acoustic oscillations using the Baryon Oscillation Spec- 
troscopic Survey data from the ninth data release of the Sloan 
Digital Sky Survey ( Ander son et al.|2012[ hereafter referred to as 
BOSS). We adopt their primary reconstructed distance constraint 
Dv(z — 0.57)/r s = 13.67 ± 0.22. Here r s is the sound horizon 
at the baryon drag epoch, and Dy (z = 0.57) is the volume element 
at a redshift z = 0.57 which depends on angular diameter distances 
and the Hubble parameter H(z). This constraint is found to be in 
excellent agreement with measurements of the distance-redshift re- 
lation from Type la supernovae (Conley et al. 2011 ; Suzuki et al. 
|2012[ >. As baryon acoustic oscillations and supernova are probing 
similar geometric properties of the Universe, using current super- 
nova data in combination with WMAP7 and BOSS yields little to 
no improvement for the majority of cosmological parameters that 
we constrain in this analysis (see Anderson et al. 2012 for more de- 
tails). We therefore do not include Type la supernovae constraints. 
We do however include a Gaussian prior on the Hubble constant, 
Ho = 73. 8±2.4kms _1 Mpc _1 , which combines constraints from 
local supernovae, Cepheid variables, and the megamaser at the cen- 
tre of NGC 4258 ( |Riess et al.|20T7] hereafter referred to as Rl 1). 



3 METHOD 

In this section we review the theory and measurement of weak lens- 
ing in tomographic redshift bins, discuss the non-linear intrinsic 
alignment model that we adopt for this analysis, and present our 
method to estimate the covariance matrix error from the H arnois-I 
|Deraps et al.| ( |2012} suite of high resolution N-body lensing simu- 
lations. We focus on a real-space shear correlation function analy- 
sis in this paper, presenting a fully 3D spherical harmonic analysis 
in |Kitching et al.H20T3) . We conclude this section describing the 
properties of the chosen tomographic redshift bins and the Popu- 
lation Monte Carlo method that we use to determine cosmological 
parameter constraints from the data. 



3.1 Weak Lensing Tomography 

Weak gravitational lensing by large-scale structure induces weak 
correlations between the observed shapes of distant galaxies. We 
parametrize galaxy shape through the complex galaxy ellipticity 
e = e\ + 1€2 - The simple relationship between ellipticity and shear, 
given in equation[T| holds for weak shear, when the ellipticity for a 
perfect ellipse with an axial ratio j3 and orientation cp, is defined as 



ei 

£2 



13+1 



1 ( cos 2<j> 
sin 2(f> 



(3) 



There are a range of different two-point statistics that have been 
proposed to extract weak lensing information from the data (see 
Schneider et al. 2002, 2010, for a comprehensive discussion of the 
relationship between these statistics). These statistics, however, all 
stem from a base measurement of the observed angular two-point 
correlation function £± which can be estimated from two redshift 
bins, i and j, from the data as follows: 



,,,, n . = E w ° w 6 bli x ^)4( x b) ± e l x (xa.)e 3 x (x b ) 

J2 w aW(, 



(4) 



Here the weighted sum, using inverse variance weights w, is taken 
over galaxy pairs with angular separation \x a — Xb\ = 0. The 
tangential and cross ellipticity parameters e t ,x are the ellipticity 
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parameters in equation ([3j rotated into the reference frame joining 
each pair of correlated objects. In this paper we only focus on this 
statistic, referring the reader to Kilbinger et al. (2013 ) who present 
a non-tomographic analysis of CFHTLenS using a wide range of 
different two-point statistics. This analysis demonstrates that for 
CFHTLenS, the cosmological parameter constraints are insensitive 
to the two-point statistic adopted for the analysis. 

The two-point shear correlation function £± (0)gg is related to 
the underlying non-linear matter power spectrum Pg that we wish 
to probe, with 



emGG = 



2n 



d££P^(£)J±(£0) : 



(5) 



where J± (£9) is the zeroth (for ) and fourth (for f _ ) order Bessel 
function of the first kind. P K (£) is the convergence power spectrum 
at angular wave number £ 



Pi j {£) = 



(6) 



where fn(w) is the angular diameter distance out to the comoving 
radial distance, w, and wh is the horizon distance. The lensing 
efficiency function, qi (w), for a redshift bin i, is given by 



q,(w) 



2c 2 a(w) 



J w 



, I , ,s!k{w' 
dw rii(w ) 



(7) 



where rii(w)dw is the effective number of galaxies in dw in red- 
shift bin i, normalized so that J rii(w)dw = 1. a(w) is the di- 
mensionless scale factor, Ho is the Hubble parameter and f2 m the 
matter density parameter at z = 0. For more details see Bartelmann 
|& Schn eider (2001 1 and references therein. 

3.2 Non-Linear Intrinsic Alignment Model 

In this paper we adopt the non-linear intrinsic alignment model 
developed by |Bridle & Kin g ( 2007 I to parametrize the contribu- 
tion of intrinsic alignments to our tomographic shear measurement. 
This model is a simplified version of the linear tidal field alignment 
model derived analytically by Hirata & Seljak (2004), based on the 
earlier work of Catelan et al. (2 001| ). |Bridle & King (2007) choose 



to make one key addition to this model by replacing the linear mat- 
ter power spectrum with a non-linear power spectrum, hence the 
name, non-linear intrinsic alignment model. This modification to 
the original model was motivated by comparisons of the model pre- 
dictions to measurements of intrinsic a lignments from data (|Man-| 
delba um"et al.|2006) , and simulations (Heymans et al. 2004J, and 
the desire to make the linear tidal field alignment model more re- 
alistic on small scales. This model has since been adopted by sev- 
eral observational analyses ( jKirk et aT. 2010; Joachim i et al.|20lT| 
Mandelbaum et al. 2 01 1| >, as it has the useful property that with 
only a single parameter A, both the II and GI contribution to the 
shear correlation function can be predicted. The non-linear intrin- 
sic alignment II and GI power spectra are related to the non-linear 
matter power spectrum as, 



Pa(k,z) =F 2 (z)P s (k,z) P GI (k,z) = F(z)P s (k,z), 



(8) 



where the redshift and cosmology-dependent modification to the 
power spectrum is given by 



F(z) = -AC lPc 



D(z) 



(9) 



Here p cr i t is the critical density at z — and D(z) is the linear 
growth factor normalized to unity today. We follow Joachimi et al. 



1 201 1 1 by parameterizing the amplitude of F{z) with a free dimen 
sionless amplitude parameter, A, and a fixed normalization con 
stent Ci = 5 x 10" 14 h~ 2 M, 



@ Mpc 3 . The value of Ci is chosen 



so that the model matches the observational results of Bro wn et al.l 
(2002) such that the fiducial model for our analysis will assume 
A = 1. In this case the GI term is negative and acts to decrease 
the overall signal. The II term, however, is always positive, inde- 
pendent of the sign of A, and acts to increase the overall signal. 
We note that F(z) differs from |Bridle & King|(2007^ , as we in- 
corporate the redshift dependent corrections to the linear tidal field 
alignment model reported in Hirata & Seljak (2010); Joachimi et al. 

poTT) . 

The II and GI contributions to the observed two-point correla- 
tion function are analogous to the GG contribution from equation|5] 



GI 



2tt 



d££CH Gl (£) J±(#), 



(10) 



with the convergence power spectrum P K replaced by the projected 
GI power spectrum Cgi or projected II power spectrum Cn, 



C u {£)- I dw y P U 



fx{v>y 



(id 



C%(£) 



qi (w KM +m(w)qj(w) 
dw P G i , , ,w 

Uk{w)] 2 \fK{w) 

(12) 

where the projection takes into account the effective number of 
galaxies n(w), and, in the case of GI correlations, the lensing effi- 
ciency qt(w) (equation [7|. Consider two non-overlapping distinct 
redshift bins, such that rii (w)rij (w) — for all w. As the II term 
comes from physically close galaxy pairs, we find for these non- 
overlapping bins Cn = 0, as expected. The GI term comes from 
the correlation of background shear with foreground intrinsic el- 
lipticities. For the same two non-overlapping bins, with the mean 
redshift of bin j greater than the mean redshift of bin i, we see 
that only the m(w)qj (w) term in the projection is non-zero, again 
as expected for GI. In practice, however, we will find that statisti- 
cal and catastrophic errors in photometric redshift estimation will 
result in some level of overlap between all the bins, so we expect 
some contribution from II and GI between all our tomographic bin 
combinations. 

Following from equation [2] we can now relate the observed 
two-point ellipticity correlation function (equation |4| to the two- 
point shear correlation function that we wish to measure (GG term, 
see equation [5} and the two types of intrinsic alignment contam- 
ination (II and GI) that we wish to marginalise over in our weak 
lensing analysis, 



(13) 



Before we continue to discuss our analysis technique to extract 
cosmological parameters from the two-point ellipticity correlation 
function, we should pause to assess how realistic the adopted non- 
linear intrinsic alignment model is. As noted in Hirata et"aL"] ( |2007| >, 
the |Catelan et al.| ( |2001| l model of linear galaxy alignment has no 
sound grounding in theory, in contrast to the analytical models of 
angular momentum correlations, and hence galaxy alignment, from 
tidal torque theory (see the review by Schafer 2009, and refer- 
ences therein). It does however provide a good match to simula- 
tions and observations, (see for example |Hirata et al.|2 007; Bridle 
& King 2007 1 and is therefore currently the favoured model. Equa- 
tion|9] assumes that the amplitude A of the intrinsic galaxy align- 
ment signal is independent of redshift, and that even at the smallest 
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physical scales, the II and GI power spectra follow the same scale- 
dependent evolution as the underlying matter power spectrum. Both 
these assumptions are clearly over-simplifications and a truly con- 
servative tomographic analysis should therefore implement a more 
flexible intrinsic alignment model, allowing for scale and redshift 
dependent perturbations around the non-linear intrinsic alignment 
model as proposed by |Bridle & King (2007). We choose not to 
take this approach, however, motivated by the null measurement of 
intrinsic alignments using a sample of galaxies that are most sim- 
ilar to the CFHTLenS galaxy population, at the median redshift 
of CFHTLenS, by Mandel baum et alT] ^201 1[ >. This observational 
analysis places an upper bound of 0.03 (2a limit) on the size of the 
systematic error induced on a measurement of as when neglect- 
ing intrinsic alignments in a CFHTLenS-like survey. We can there- 
fore conclude that the systematic error introduced by the presence 
of intrinsic alignments in CFHTLenS is likely to be sufficiently 
small that any additional systematic error we introduce by using 
an inflexible single-parameter intrinsic alignment model will be in- 
significant compared to our statistical errors (which we show in 
Section |4~2"l result in a la error of ±0.04 on as at fixed Q m , when 
intrinsic alignments are marginalised over). This rationale will no 
longer hold for future surveys, however, where the large survey ar- 
eas will reduce the statistical errors, highlighting the pressing need 
for improved models and observations of intrinsic galaxy align- 
ments. 



3.3 Covariance Matrix Estimation 

One of the major challenges with a tomographic lensing analysis of 
the type that we use in this analysis, is the construction of a suffi- 
ciently accurate covariance matrix that accounts for the significant 
levels of correlation between both the angular and redshift binned 
data. For N t tomographic redshift bins, Ng angular scales and con- 
sidering both the £^ and ^ components of the shear correlation 
function between redshift bins i and j, we have a total number of 
data points p given by 

p = NeN t (N t + l). (14) 

It is both the data vector, D{p), and the inverse of the correspond- 
ing p x p covariance matrix, C ~ , that we use in the cosmological 
parameter likelihood analysis. In this section we show that it is the 
covariance matrix estimation that limits the maximum number of 
data points that we can accurately analyse, motivating and justify- 
ing our choice of p — 210 in Section [3~4| 

There are a number of methods that can be used to estimate a 
covariance matrix. Bootstrap and Jackknife techniques can be used 
to estimate the data covariance directly from the survey ( |Wall &| 
|Jenkins|2003) , and various de-noising techniques can be applied to 
reduce the impact of noise bias when the matrix is inverted (see for 
example Norberg et al. 2009). Analytical functions can be derived 
for random Gaussian convergence fields { [Schneider et~aLl|2"002"l 
Kilbinger & Schneider 2004) and scale-dependent non-Gaussian 
corrections can be applied to those analytical functions as deter- 
mined through comparison to N-body lensing simulations ( |Sem-| 
boloni et al. 2007, Sato et al. 2011 1. Analytical functions can also 
be derived for log-normal convergence fields or using a halo model, 
which are found to be a significantly better approximation com- 
pared to a Gaussian field. These analytical functions, however, still 
require calibration using N-body simulations (Hilber t et al.|20lT| 
|Kayo et al.[ 2012). With sufficient numbers of 3D N-body lensing 
simulations, the covariance can also be estimated directly from the 



simulations. These can be populated to emulate the survey geome- 
try, masks, redshift distribution, galaxy number density and intrin- 



sic ellipticity and shape measurement noise for th e survey {H ubert 



etal. 


2009 


|Vafaei et al.|2010; Kiessling et al.|201 1 ||Harnois-Deraps 


etal. 


2012 


i. 



There are advantages and disadvantages to each method. Di- 
rect estimation of the covariance matrix from the survey provides 
a cosmology-independent, but noisy estimate. In the inversion, the 
covariance matrix can therefore become unstable. Furthermore, this 
method is only suitable when the survey contains many indepen- 
dent lines of sight, adequate to evaluate sampling variance errors 
which dominate at large scales. As such, this direct estimation is 
only suitable for very large area surveys. Analytical functions are 
precise, but approximate. They can be calculated for any cosmol- 
ogy, but require cosmology and scale-dependent corrections, for 
angular scales < 20 arcmin, calibrated with N-body simulations, 
to account for the Gaussian or log-normal approximations made in 
their calculation. N-body simulations are costly, and for the real- 
space correlation analysis presented in this paper, the requirements 
on simulation resolution and size is demanding. Low particle res- 
olution in the simulation results in an artificial lack of power that 
propagates to all scales with real-space statistics. The finite simu- 
lation box size also truncates density perturbations on scales larger 
than the box length which results in a suppression of power in the 
large-scale real-space dark matter halo correlation function mea- 
sured from the simulations (Power & Knebe 2006). High resolu- 
tion, large area, 3D lensing simulations that are suitable for our 
real-space analysis are therefore rare, and for a fixed cosmology 
( |Hilbert et al.|2009||Harnois-Deraps et al.|2012| >. 

In this analysis we estimate a covariance matrix from 
the three-dimensional N-body numerical lensing simulations of 
|Harno is-Deraps et al. (2012). Light cones are formed from line of 
sight integration through independent dark matter particle simula- 
tions. The simulated cosmology matches the 5-year WMAP flat 
ACDM cosmology constraints from Dun kley et al.| ( |2009| >. The 
1024 3 particle simulations have a box size of 147.0h- 1 Mpc or 
231.1 h _1 Mpc, depending on the redshift of the simulation. The 
boxes are grafted in projection such that each high resolution line 
of sight simulation has a real space resolution of 0.2 arcmin in the 
shear field, spanning 12.84 square degrees sampled at 26 redshift 
slices between < z < 3. The two-point shear statistics measured 
in real-space from the simulations closely matches the theoretical 
predictions of the input cosmology from 0.5 < 9 < 40 arcmin 
scales at all redshifts (Harnois-Deraps et al. 2012). We therefore 
limit our tomographic analysis to these angular scales where we 
can obtain an accurate estimate of the covariance matrix. 

We populate each simulation by mapping 12.84 square degree 
sections of the galaxy spatial distribution from the survey on to 
the simulated shear field, and the corresponding redshift distribu- 
tion and galaxy weights. With this method the survey masks and 
effective source density are correctly replicated in the simulations. 
We create a redshift realization of each line of sight, by sampling 
a galaxy redshift z at random from each galaxy P(z) in the sur- 
vey. We verified that our results were insensitive to this realiza- 
tion method by repeating the measurement of the covariance for 
different random seeds. A galaxy shear 7 is then assigned at that 
redshift by linearly interpolating between the shear output at dis- 
crete redshifts by the lensing simulations. Intrinsic ellipticity and 
shape measurement noise is added by randomly sampling from a 
zero-mean Gaussian of width a e = 0.279 per ellipticity compo- 
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nent as measured from the dataj^] Note that we found no signifi- 
cant variation in the shot noise <r e when measured as a function 
of photometric redshift within our high confidence redshift range 
0.2 < zbpz < 1.3. This was determined by binning the data into 
photometric redshift slices of width Azbpz = 0.1, where a e for 
each bin was found to lie within two standard deviations of the 
mean of the full sample. 

|Harnois-Deraps et al.| ( |2012"| > create a total of 184 fully inde- 
pendent N-body simulations. In order to gain a sufficient number of 
realizations such that the correction for correlated noise in the in- 
verse covariance matrix will not bias our results (see the discussion 
in Section [3.3.1| l, we are therefore required to split each of the sim- 
ulations into 9 semi-independent simulations. Each simulation then 
spans 1 .4 square degrees such that the centres are separated by more 
than twice the largest angular scale probed by this analysis. This 
method to increase the number of simulations by sub-division, also 
employed by Schrabb ack et al.| ( |2010| l, yields n M — 1656 semi- 
independent line of sight simulations. In Section [4~2| we show that 
the low level of correlation expected between each group of 9 sim- 
ulations, introduced by splitting each fully independent simulation 
into sub-fields, does not impact significantly on our results. The 
data vector Z? M of the p components of (8), is then measured 
from each simulation line of sight ra M and the covariance matrix is 
estimated, 



C(M) 



A s n, 



^=1 



D a )(D£-D b 



(15) 



where D is the average of the data measured over the total n M 
simulations. The area scaling term A s accounts for the difference in 
area between one line of sight in the simulation and the CFHTLenS 
survey area used in the analysis I Schn eider et al.|2 002 1. 

The covariance matrix we derive with this method ignores the 
impact of intrinsic galaxy alignments. As our fiducial A = 1 intrin- 
sic alignment model is expected to lower the amplitude of the ob- 
served two-point correlation function only at the few per cent level, 
a GG-only covariance matrix is therefore expected to only slightly 
over-estimate our errors^] We therefore conclude a GG-only co- 
variance is sufficient for our analysis supported by the findings of 
|GrocuttH2012) who analysed the covariance matrix derived from 
simulated II+GI+GG Gaussian fields created using the method of 
|Brown &"Tj attye (2011). For future tomographic analyses we will 
populate the N-body lensing simulations used in this analysis with 



2 |Miller et al.||20lj) show that the true intrinsic galaxy ellipticity distribu- 
tion is not well approximated by a Gaussian distribution. For the purposes 
of evaluating the dominant shot noise component of the covariance matrix, 
however, we argue that a Gaussian approximation is sufficient. The width 
of the Gaussian model used is calculated directly from the data such that the 
variance of the data and simulated ellipticity distributions are the same, even 
if the shapes of the distributions differ. For the sub-dominant mix-term part 
of the covariance, however, which arises from correlations between galaxy 
ellipticities and cosmic shear, our Gaussian approximation could lead to 
a mild underestimate of the covariance and this should be investigated in 
future work. 

3 In the analysis that follows we find A to be consistent with zero for 
the full galaxy sample and a sample of late-type galaxies, supporting this 
strategy. For early-type galaxies we find a significant value for A. In this 
case, however, shot noise dominates the covariance matrix as the early-type 
galaxy sample comprises only 20 per cent of the full galaxy population. In 
this case it is therefore again a reasonable strategy to ignore the impact of 
intrinsic galaxy alignments on the covariance matrix. 



an intrinsic ellipticity component (see for example Hey ma ris et al.| 
|2006l >. 

The covariance matrix we derive with this method is depen- 
dent on the fixed simulation cosmology. |Eifler et al.| (|2009) inves- 
tigate the impact of including a varying cosmology covariance ma- 
trix in a cosmic shear analysis, concluding that fixing the cosmol- 
ogy of the covariance matrix has a non-negligible impact on the 
size of the likelihood contours. This is an effect that becomes less 
pronounced as the statistical power of a survey increases. Kilbinger 
et al. (2013 ) therefore present a detailed analysis of the impact of 
including a varying cosmology covariance matrix in the analysis of 
the 2D cosmic shear measurement for the CFHTLenS data. They 
also present a novel method to combine the advantages of all the 
covariance matrix estimation methods to create a cosmology de- 
pendent covariance matrix. On scales 6 < 30 arcmins, the non- 
Gaussian sampling variance is estimated from the N-body lensing 
simulations of Harnois-Deraps et al. ( 2012 ), as described above, but 
setting a e = 0. A varying cosmology estimate is then obtained by 
scaling the data vector Z? M by the ratio of the desired cosmology 
to the fixed simulation cosmology. For the larger angular scales, 
where the finite box size causes power to be underestimated in the 
simulations, a Gaussian analytic model is used, which has been 
shown to be a good approximation on these scales ( Sembolo ni et al.| 
2007; Sato et al. 201 1 ). The shot-noise term of the covariance ma- 
trix is also estimated analytically (Schneider et al. 2002) and the 
cosmology variation in the mixed term, which accounts for the co- 
variance between the cosmic shear and shot noise, is modeled using 
the |Eifler et al.| (|2009) fitting formulae. The conclusion of this de- 
tailed covariance matrix analysis of Kilbinger et al. (2013) is that, 
for a CFHTLenS-like survey, the impact of using a fixed cosmol- 
ogy covariance matrix, or a varying cosmology covariance matrix 
is marginal, particularly when CFHTLenS results are used in com- 
bination with other cosmology surveys. Kilbi nger et aTJpOB) also 
show a good agreement between the diagonal components of the 
covariance matrix with a lackknife estimate of the covariance from 
the data. We therefore conclude that using the Harnois-Deraps et al. 
l |2012} simulations to produce a fixed cosmology covariance ma- 
trix, as estimated using Equation[l5j is sufficiently accurate for our 
purposes. 

3.3.1 Inverse Covariance Estimation 

For the likelihood analysis of the CFHTLenS data, we require the 
inverse of the covariance matrix. Whilst we consider our measure- 
ment of C from N-body lensing simulations to be an unbiased es- 
timator of the true covariance matrix C, it will have an associated 
measurement noise from averaging over a finite number of semi- 
independent realizations n M . This measurement noise means that 

C is not an unbiased estimate of the true inverse covariance ma- 
trix ([Hartlap et al. 2007 ). Assuming Gaussian measurement errors 

on C, an unbiase d estimate of the true inverse covariance matrix is 

-l 



derived in 



Anderson 



(2003 



QA = 



where C 

-P 



CtA I 



and 



(16) 



Hartlap et al. 1 2007 1 shows that for p/n M < 0.8, this correction 
produces an unbiased estimate of the inverse covariance matrix 
C 1 and we use this correction in our analysis. 

Unfortunately, even though the use of the |Anderson| (2003] ) 
correction creates an inverse covariance that is now unbiased in the 
mean, there are still associated measurement errors on the individ- 
ual components of C which now become boosted by a factor a a ■ 



© 0000 RAS, MNRAS 000, 000-000 



8 C. Heymans et al. 



Hartlap et al. ( 2007]) show that boosting this measurement noise in 
the inverse covariance, impacts on the size of the confidence re- 
gions resulting from a likelihood analysis. They find that the area 
of the confidence regions can erroneously grow by up to 30 per 
cent as p/n M — > 0.8. The consequence is that for high values of 
p/rifi, parameter constraints appear less significant than they truly 
are. If we require that the inclusion of the Anderson ( 2003 1 correc- 
tion does not increase the area of our Bayesian confidence regions 
by more than ~ 5 per ce nt, |Hartlap et al7| |2007 1 estimate a limit 
of p/n^ < 0.12 (see also |Grocutt|20 1 2 for similar conclusions). 
We note that the tomographic analysis of the COSMOS survey by 
|Schrabback et al.| ( |2010} estimated the covariance matrix from the 
Millennium Simulation with p/n M = 0.55. Whilst well within the 
requirements for an unbiased Anderson (2003 ) correction, Hartlap 
et al. (2007) predicts such a correction would increase the area of 
the confidence regions in a likelihood analysis by ~ 12 per cent. 
As such the Schrabback et al. (2010) cosmological parameter con- 
straint errors are expected to be slightly overestimated. 

As we have a fixed number of N-body lensing simulations 
which determines it is the accuracy that we require for the 
inverse covariance matrix C 1 that sets the maximum number of 
data points p in our analysis. The number of tomographic bins N t 
and angular scales N# is therefore set by the number of N-body 
simulations that we have at our disposal. For p/n M < 0.12, and 
n M = 1656, (see Section [3~3| , we should therefore limit our analy- 
sis to p < 200. 



3.4 Tomographic analysis and redshift distributions 

In a tomographic weak lensing analysis there is always a choice 
to be made for the number of tomographic redshift bins, N t , and 
the number of scales probed, in our case angular scales, Na. As 
the number of redshift and angular bins is increased, the amount of 
information increases. A saturation limit is eventually reached be- 
yond which the data points become so correlated that the extra in- 
formation gained with each incremental increase in the number of 
bins becomes marginal. With an unlimited number of N-body lens- 
ing simulations from which to make an unbiased covariance matrix 
estimate, the optimal number of tomographic bins will depend on 
the photometric redshift accuracy of the survey, and the method by 
which the contamination from intrinsic galaxy alignments is miti- 
gated in the analysis. Bridle & King (2007) show that for a survey 
with a photometric redshift scatter of a z = 0.05(1 + z), using 
N t ~ 8 brings the cosmological parameter constraints to within 20 
per cent of the best attainable with a fully 3D approach. This is in 
contrast to the conclusions of earlier cosmic-shear only optimiza- 
tions, which found Nt ~ 3 to be optimal ( Simon et al. 2004) |Ma| 
|et al.|20 06 l. This difference indicates the importance of using finely 
binned tomographic redshift slices when mitigating intrinsic align- 
ment effects. [Grocutt| f2012 1 investigate the dependence of cosmo- 
logical parameter constraints when varying the number of tomo- 
graphic redshift bins, N t , and the number of angular scales probed, 
Ne, simultaneously. A non-linear intrinsic alignment model was 
assumed for the II and GI contamination (see Section [3~2] l. In this 
analysis the cosmological parameter constraints were found to be 
less sensitive to increases in Ne, in comparison to increases in 
Nt. This is expected for the single-parameter non-linear intrinsic 
alignment model, as the cosmic shear, GG, and non-linear intrinsic 
alignment II and GI power spectrum, vary smoothly with scale and 
the relative amplitude between the II, GI and GG power for each 
redshift bin is fixed as a function of scale. As the number of data 



Table 1, Tomographic redshift bin selection. Galaxies are selected based on 
their maximum posterior photometric redshift estimate zbpz- The median 
redshift z m and mean redshift z for each bin is calculated from the effective 
redshift distribution as measured by the weighted sum of the photometric 
redshift error distributions P(z). 



Bin 


2 BPZ 




z 


1 


0.20 - 0.39 


0.28 


0.36 


2 


0.39 - 0.58 


0.48 


0.50 


3 


0.58 - 0.72 


0.62 


0.68 


4 


0.72 - 0.86 


0.82 


0.87 


5 


0.86 - 1.02 


0.93 


1.00 


6 


1.02 - 1.30 


1.12 


1.16 



points p scales as N t (Nt + 1), however, even small increases in Nt 
can quickly lead to an unstable covariance matrix. 

Motivated by the findings of Brid le & King| ( |2007| > and |Gro-| 
|cutt| ( 12012"} , and with the limitation that the total number of data 

. we choose to use N± 



points p < 200 (see Section 3.3.1 1, we choose to use N t = 6 red- 
shift bins and Ne = 5 angular bins such that our total number of 
data points p = 210. The angular range is chosen to be spaced 
equally in log(#) between 1.5 < 9 < 35 arcmin, where the maxi- 
mum angular scale is determined by the limitations of the N-body 
lensing simulations used to determine the covariance matrix. We 
select the Nt = 6 redshift bins to span our high confidence redshift 
range 0.2 < zbpz < 1-3 such that the effective angular num- 
ber density of galaxies in each redshift bin is equal. The effective 
number density includes the shear measurement weights w such 
that the intrinsic ellipticity noise in each bin is equal. This choice 
is in contrast to a cosmic shear signal-to-noise optimised redshift 
bin selection which would lead to much broader bins at low red- 
shift. Such optimization is undesirable for our purposes, as it is 
the lowest redshift bins where the presence of intrinsic alignments 
is most prominent. Table[T|lists the resulting redshift selection for 
each tomographic bin. The median redshift z m and mean redshift 
z is calculated from the effective redshift distribution as measured 
by the weighted sum of the photometric error distributions P(z). 
These error distributions extend out to zbpz = 3.5 which skews 
the mean redshift measurement, relative to the median, particularly 
in the lowest redshift bin. 

Figure [T] compares the effective redshift distribution for each 
tomographic bin as determined from the maximum posterior red- 
shift zbpz (upper panel) and by the weighted sum of the photomet- 
ric redshift error distributions P(z) (lower panel). In order to high- 
light the differences between the redshift distributions measured 
from these two different estimators, the binning in the upper panel 
is chosen to be significantly finer than the typical CFHTLenS pho- 
tometric redshift error a z ~ 0.04(1 + z) jHildebrandt et al.|2012| >. 
The fine structure revealed by this binning therefore illustrates red- 
shift focusing effects arising from the photometric redshift mea- 
surement, not true physical structures. Accurate measurements of 
P(z) for each galaxy allows us to fully account for these focussing 
effects, in addition to overlapping redshift distributions and catas- 
trophic redshift outliers in our analysis (see Benjamin et al. 2013 
for detailed analysis of the P(z) used in this analysis). For our in- 
trinsic alignment analysis it is particularly important to quantify the 
degree of overlap between redshift bins as the II term is only signif- 
icant for physically close galaxy pairs. It is therefore the summed 
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Figure 1. Tomographic redshift distribution. The upper panel shows the 
effective weighted number of galaxies as a function of their maximum pos- 
terior photometric redshift estimate, separated into six tomographic bins 
between 0.2 < zbpz < 1-3. The effective weighted number of galaxies in 
each redshift bin is constant. The lower panel shows the redshift distribution 
for each selected bin as estimated from the weighted sum of the photometric 
redshift probability distributions P(z). 



P(z) redshift distributions displayed in the lower panel of Figure[T] 
that we use in this analysis. 



3.5 Population Monte Carlo Sampling likelihood analysis 
method 

In this study we perform a Bayesian likelihood analysis of 
CFHTLenS and the auxiliary data, discussed in Section [2] to con- 
strain the parameters of a range of cosmological models. To calcu- 
late the likelihood values we u se the Population Mont e Carlo sam- 
pling software CQSMOPMCj^lKiibin ger et al.|201l| , modified to 
include an optional simultaneous fit of cosmic shear and the intrin- 
sic alignment model outlined in Section [3~2| Future releases of this 
software package will include this option. The Population Monte 
Carlo method is described in Wraith et al. ( 2009 1 along with a com- 
parison to the more standard Markov-Chain Monte Carlo method 
for cosmological parameter estimation. We also refer the reader to a 
detailed discussion of the COSMOPMC analysis of 2D CFHTLenS 
cosmic shear data in Kilbinger et al. (2013 1 for further information 
about the methodology. 

We assume a matter power spectrum derived fro m the | Eisen-| 
|stein & Hu|ll998| > transfer function with a |Smith et aT1 ( |2003) norT 
linear correction. For dark energy cosmologies, where the equa- 
tion of state of dark energy parameter, wo 7^ —1, a modulation 
of the non-linear power is required ( McDonald et al. 2006 1 which 
we apply using of the scaling correction from [Schrabback et al.| 



p0T0| ; Ref regier etall pOTT) . The |Smith et"aL] ( |2003| > halo-model 
prescription for the non-linear correction has been calibrated on 
numerical simulations and shown to be accurate to between 5 per 
cent and 10 per cent over a wide range of k scales ( |Eifler|201 1) 
and found to be of sufficient accuracy for the statistical power of 
CFHTLenS (Vanderveld et al. 20121. Whilst our assumed transfer 
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function includes baryonic oscillations on large scales, we are un- 
able to include the uncertain effects of baryons on small physical 
scales. Semboloni et al7|l |201 1| present an analysis of cosmological 
hydrodynamic simulations to quantify the effect of baryon physics 
on the weak gravitational lensing shear signal, using a range of dif- 
ferent baryonic feedback models. For the angular scales we use 
we would expect baryons to induce at most a ~ 10 per cent de- 
crease in the signal relative to a dark matter only Universe, in the 
mid-to-high redshift tomographic bins where our highest signal-to- 
noise measurements are made. This is assuming the 'AGN feed- 
back' model which leads to the largest changes in the matter power 
spectrum of the simulations that were studied by Semboloni et al. 
l |2011) , where we note that this scenario is the one that matches 
observed gas fractions in groups. In the cosmological analysis that 
follows, we present an additional conservative analysis where the 
tomographic data most susceptible to significant errors caused by 
baryonic or non-linear effects are removed (see Benjami n et al.| 
|2013| for further discussion). If significant errors exist, however, 
the inclusion and marginalisation over the intrinsic alignment am- 
plitude A in our analysis, which modulates the amplitude of the 
observed shear power spectrum, should work to some extent, to 
reduce the impact of these effects in addition to mitigating contam- 
ination by intrinsic galaxy alignments. 

We use COSMOPMC to analyse CFHTLenS and WMAP7 in- 
dependently. For the combined results with BOSS and our assumed 
Ho prior from Rl 1, we importance-sample the WMAP7-only like- 
lihood chain, multiplying each sample point with the CFHTLenS, 
BOSS and Rll posterior probability. For our CFHTLenS-only flat 
ACDM analysis we limit our parameter set to the matter density 
parameter, f2 m , the amplitude of the matter power spectrum con- 
trolled by erg, the baryon density parameter tli,, the Hubble param- 
eter h, and the power spectrum spectral index n s . With WMAP7 we 
also include into the parameter set the reionisation optical depth r, 
the Sunyaev-Zel'dovich template amplitude Asz, and the primor- 
dial amplitude of the matter perturbations A| , from which we de- 
rive erg. The equation of state of dark energy parameter, wo and 
dark energy density parameter Q^c are also included for non-flat or 
non-A cosmological models. We use flat priors throughout which 
are broad enough to cover the full 3<j posterior distribution in each 
parameter direction for each combination of data. Throughout the 
paper we quote and plot 68 per cent and 95 per cent Bayesian con- 
fidence or credibility regions. These regions contain 68 per cent 
and 95 per cent of the posterior probability determined from the 
multi-dimensional distribution of points from the PMC parameter 
sample. All figures showing the resulting joint-constraints on two 
parameters, are marginalised over the multi-dimensional parameter 
space that is not shown. 



4 RESULTS 

Figure [2] presents the observed two-point correlation function 
f V (9) for every tomographic bin combination in our chosen six 
redshift bin analysis. With N t tomographic bins, there are N t (Nt + 
l)/2 independent combinations, or 21 combinations in our case. 
The panels show the different ij bin combinations, ordered with 
increasing redshift bin i from left to right, and increasing redshift 
bin j from lower to upper, where the redshift distributions of each 
bin are shown and tabulated in Section [3~4l The auto-correlated 
bins lie along the diagonal. The data points are calculated using 
the shear correlation function estimator in Equation [4] correlat- 
ing pairs of galaxies within the full mosaic catalogue for each of 
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Figure 2. The observed two-point correlation function £_P (8). The panels show the different ij redshift bin combinations, ordered with increasing redshift bin 
i from left to right, and increasing redshift bin j from lower to upper. Refer to table[T]for the redshift ranges of each tomographic bin. The errors are estimated 
from an analysis of N-body lensing simulations as discussed in Section [33] The theoretical curves show our fiducial total GG+GI+II signal as a solid line. 
When distinguishable from the total, the GG only signal is shown dashed. The magnitude of the GI signal is shown dot-dashed (our fiducial GI model has 
a negative anti-correlated signal) and the II signal is shown dotted, where the amplitude is more than 10 — 7 . The results of the broad two-bin tomographic 
analysis of Benjamin et al. 1 2013 1 are shown in the lower right corner. 



the four CFHTLS fields. The measurements from each field are 
then combined using a weighted average, where the field weight is 
given by the effective number of galaxy pairs in each angular bin. 
Note that the results for each ij bin from each field were found 
to be noisy but consistent (see Kilbi nger et al. 12013] for measure- 
ments of the higher signal-to-noise 2D shear correlation function 
for each CFHTLS field). The errors, which include sample vari- 
ance, are estimated from an analysis of N-body lensing simulations 
as discussed in Section [3~3l We remind the reader that the data are 
highly correlated, particularly in the low redshift bins. The theoreti- 
cal curves show our fiducial WMAP7 best-fit cosmological param- 
eter model, with an A — 1 non-linear intrinsic alignment model, to 
be a good fit to the data. A possible exception to this are data from 
tomographic bin combinations that include the lowest redshift bin, 
which we discuss further in Section pTT] The individual components 
are shown; GG (dashed), GI (dot-dashed) and II (dotted) models 
with the total GG+GI+II shown as a solid line. For comparison we 
also show the results of the broad two-bin tomographic analysis of 



|Benjamin et al.| ([2013 1 in the lower right corner to demonstrate the 
low-level of II and GI contamination expected for this high redshift 
selected analysis. 



4.1 Tomographic Data Visualization 

With 21 tomographic bin combinations, two statistics £1(0) and 
C-(0), and 5 angular scales, we have a total of p = 210 data 
points, half of which are shown in Figure [2] In the cosmological 
parameter constraints that follow, it is this large data vector, and a 
correspondingly large covariance matrix, that we use in the likeli- 
hood analysis. Purely for improving the visualization of this large 
data set, however, we propose the following method to compress 
the data, motivated by the different methods of Massey et al.|(2007| > 
and |Schrabback et al.|p010) . 

To compress angular scales, we first calculate a WMAP7 cos- 
mology GG-only theory model t^ d for each redshift bin combina- 
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tion ij and each statistic (+/—). We then define a free parameter 
a± which allows the overall amplitude of the model to vary, but 
keeps the angular dependence fixed. The best-fitting amplitude a± 
is then found from a \ 2 minimization of a^^ d (6) to the shear 
correlation functions measured at 5 angular scales in each ij bin 
and each statistic. A best-fitting value of = 1 implies the data 
in bin ij are well-fit by a WMAP7 GG-only cosmology. Following 
|Schrab back et al. (2010), each bin is then assigned a single value 
of a 13 Cid(® = 1) which can be interpreted as the amplitude of 
the two-point shear correlation function measured in bin ij at an 
angular scale of 8 = 1 arcmin. 



To compress the information in the redshift bin combination, 
we calculate the lensing efficiency function qi(w) (equation |7J for 
each redshift bin i, and then determine the peak redshift z pea k of 
the combined lensing sensitivity qi(w)qj(w) for each redshift bin 
ij combination. This peak redshift locates the epoch that is the 
most efficient at lensing the two galaxy samples in the redshift 
bin combination ij, but we note that these distributions are very 
broad, particularly for the redshift bins with a significant fraction 
of catastrophic outliers in the photometric redshift distribution (see 
Figure [TJ. 



Figure [3] shows the resulting compressed 21 data points for 
each statistic, (circles) and (crosses), plotting a tJ £g J d (# = 
1') against z poa k. This can be compared to the fiducial cosmology 
prediction (shown dotted, by setting a = 1). Note that the rela- 
tively high fraction of catastrophic redshift outliers in the lowest 
redshift bin impacts on the expected signal measured from redshift 
bin combinations that include this bin. The expected increase in 
signal, as 2 pea k increases, is therefore not smooth. This can be seen 
in the theoretical curve in Figure [3] which displays a slight kink at 
-Zpcak — 0.22. To recover a lJ from this figure, one simply divides 
the value of each data point by the value of the fiducial model, 
shown dotted, at that z pC ak. Consistent values for a lJ are measured 
from both the £ + and £_ statistic. We find a signal that rises as 
the peak redshift of the lensing efficiency function increases; the 
more large scale structure the light from our background galaxies 
propagates through, the stronger the lensing effect. In general, the 
data are well-fit by the WMAP7 GG-only fiducial model, but we 
do see an indication of an excess signal at low redshifts where, 
for a fixed angular scale, the smaller physical scales probed are 
more likely to be contaminated by the intrinsic galaxy alignment 
signal. This is however also the regime where the analysis is most 
affected by catastrophic outliers in our photometric redshift distri- 
bution. Based on the cross-correlation analysis of Benjamin et al. 
( |20 13| > we expect these errors to be accounted for by our use of 
photometric redshift distributions P(z). In Heymans et al.|p0l2| , 
we also show that the catalogues used in this analysis present no 
significant redshift-dependent systematic bias when tested with a 
cosmology-insensitive galaxy-galaxy lensing analysis. This gives 
us confidence in the robustness of our results at all redshifts. We 
note that in order to make this visualization of the data, the differ- 
ent redshift bin combinations and the and statistics are con- 
sidered to be uncorrelated. The plotted la errors on a are therefore 
underestimated but we re-iterate at this point that this data com- 
pression is purely for visualization purposes and it is not used in 
any of the cosmological parameter constraints that follow. 




UfT-fl 



_i_ 



0.2 0.3 0.4 0.5 0.6 
z peak 

Figure 3. Compressed CFHTLenS tomographic data where each point rep- 
resents a different tomographic bin combination ij as indicated by z pea k> 
the peak redshift of the lensing efficiency for that bin combination. The 
best-fitting amplitude a' J of the data relative to a fixed fiducial GG-only 
cosmology model is shown, multiplied by the fiducial model at 6 = 1 
arcmin for £4. (circles) and §_ (crosses, offset along the z pca k ax is f° r 
clarity). The error bars show the la constraints on the fit. The data can be 
compared to the fiducial GG-only model, shown dotted. Note that the colour 
of the points follow the same colour-scheme as Figure[TJ and indicates the 
lower redshift bin that is used for each point. 



4.2 Comparison of parameter constraints from weak lensing 
in a flat ACDM cosmology 

The measurement of cosmological weak lensing alone is most sen- 
sitive to the overall amplitude of the matter power spectrum. This 
depends on a degenerate combination of the clustering amplitude 
as and the matter density parameter f2 m , and it is therefore in this 
parameter space that we choose to compare the constraints we find 
from weak lensing alone using different analysis techniques. We 
limit this comparison to flat ACDM cosmologies. Figure [4] com- 
pares three cases. In blue we show the 68 per cent Bayesian confi- 
dence limits from a 2D weak lensing analysis of CFHTLenS, lim- 
ited to the same angular scales as our tomography analysis with 
9 < 35 arcmin. This can be compared to the 68 per cent con- 
straints from our 6-bin £± tomographic lensing measurement when 
intrinsic alignments are assumed to be zero (pale blue) and when 
the amplitude of the intrinsic alignment model is allowed to be a 
free parameter and is marginalised over (pink). All three measure- 
ments are consistent and can be compared to the best-fit WMAP7 
results shown as a black cross for reference. 

Table [2] lists the parameter constraints, for the three cases 
shown in Figure|4] on the combination crg(S7 m /0.27) a . The param- 
eter a is derived from a fit to the likelihood surface to determine the 
direction that is orthogonal to the ag — f2 m degeneracy direction. 
These results can be compared to the 2D CFHTLenS constraints 
from Kilbinger et al. (20131, where large angular scales were in- 
cluded in the analysis, and a 2-bin tomography analysis from |BerT| 
[jamin et aL|^2013| , limited to the same angular scales considered in 
this analysis. We find excellent agreement between the cosmolog- 
ical results from the different analyses, indicating that ignoring in- 
trinsic alignment contamination in Kilbinger et al. ( 2013 1 and Ben- 
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Figure 4. Flat ACDM parameter constraints (68 per cent confidence) on the 
amplitude of the matter power spectrum controlled by erg and the matter 
density parameter fi m from CFHTLenS-only, comparing three cases: 2D 
weak lensing (blue) and 6-bin tomographic lensing where intrinsic align- 
ments are assumed to be zero (pale blue) and are marginalised over (pink). 
For reference, the black cross shows the corresponding best-fit values from 
WMAP7 jKomatsuetal.|2011l 



Table 2. Constraints orthogonal to the erg ~ Q m degeneracy direction for a 
number of different types of lensing analyses: 2D weak lensing to 8 < 350 
aremm (Kilbinger et al.|2013> and to 8 < 35 arcmin, 2-bin tomographic 
lensing I Benjamin et al. 2013 1, and 6-bin tomographic lensing where intrin- 
sic alignments are assumed to be zero or are marginalised over (our primary 
result indicated in bold). We also present constraints for two conservative 
analyses to test the covariance matrix and our sensitivity to potential error 
in the assumed non-linear correction to the matter power spectrum (see text 
for more detail). 



Data 



<7 S (n m /0.27) Q 



2D Lensing 
8 < 350' 
8 < 35' 



Kilbinger et al. 



2013 



0.59 ±0.02 
0.61 ± 0.02 



0.787 ±0.032 



0.791 



+0.052 



2-bin tomography: 



< 40' Benjamin et al 



2013 



0.55 ±0.02 



0.771 ± 0.040 



6-bin tomography (all 8 < 35'): 
A = 

A marginalised 

n-fi = 736 covariance matrix 
Low 8 scales removed 



0.52 ±0.02 

0.46 ± 0.02 0.774 

0.48 ± 0.03 0.768 

0.45 ± 0.03 0.774 



fl 7RQ+0.024 
u.(OO_ 032 



+ 0.032 
-0.041 

0.032 
0.041 
0.038 
0.057 



[jamin et al. (20131 did not introduce any significant bias in their 
results. Differences in the values of the a parameter arise from a 
number of factors. The strength of the lensing signal is modulated 
by Q m in a manner which is sensitive to both the source redshift 
distribution and the angular scales under consideration ( Bernardeau 
|et al.|1997) . Furthermore the degeneracy contours are not perfectly 
represented by a power law, so the value of a is not our key interest 
here. 

Focusing first on the constraints from tomography and 2D 
lensing limited to the same angular scales but ignoring intrinsic 
alignments (shown blue and pale blue in Figure]?), we find close to 
a factor of two improvement in the constraint on <Tg(S7 m /0.27) a , in 
addition to an improvement in degeneracy breaking between erg and 
S7 m , when tomographic bins are considered. Unfortunately, how- 
ever, our tomographic analysis is limited by the extent of the N- 
body simulations used to determine our covariance matrix, which 
forces us to lose the large angular scales considered in the 2D anal- 
ysis from Kilbin ger et al.|j2013^ . Comparing the constraints from 
tomography limited to 9 < 35 arcmin, with 2D lensing out to 
6 = 350 arcmins, we find similar constraints on erg(f2 m /0.27) a . 
This demonstrates how, in this parameter space, the large angular 
scales are adding as much information in a 2D Lensing analysis 
as the additional redshift bins add in a tomographic analysis of the 
same data. This motivates future work to remove the current an- 
gular limitations imposed by the tomographic covariance matrix 
estimation method that we use in this analysis. 

Adding in the intrinsic alignment model, and hence an ad- 
ditional free parameter A, broadens the parameter constraints, as 
expected, reducing the constraining power on as(Qm/0-27) a by 
roughly 30 per cent compared to a GG-only analysis. For fl m = 
0.27, however, we find very little difference in the best-fit value 
of erg which changes by 0.01. Larger deviations between the two 
analyses are however seen for higher values of O m . Figure|4]shows 



that the 68 per cent confidence region shifts to slightly lower Sl m 
and higher erg. For a fixed fi m , the resulting best-fit erg is there- 
fore slightly lower when intrinsic alignments are marginalised over. 
This behaviour is unexpected for a conventional intrinsic alignment 
model where the negative GI signal dominates over the positive II 
signal such that the total GG+GI+II signal observed is less than the 
GG signal alone. For the fiducial A — 1 intrinsic alignment model, 
a GG-only analysis would therefore underestimate erg for a fixed 
f2 m . For this CFHTLenS-only analysis, however, we instead find 
a preference for a negative value of A — —1.601} g4, and hence 
the GG-only analysis prefers higher values of ag . We explore and 
discuss this result in more detail in Section|5]but re-iterate that the 
differences we have commented upon here are well within our 2cr 
errors and are therefore not significant. 

Finally we perform two conservative analyses to ensure the 
robustness of our results, the constraints from which are reported 
in the lower two rows of table[2]. The first is to compare constraints 
when we use a covariance matrix constructed from n M = 736 
semi-independent lines of sight (where each fully independent N- 
body lensing simulation is split into 4 sub-simulations) instead of 
the standard 7i M = 1656 analysis that we use throughout this pa- 
per. The excellent agreement between the results from the two esti- 
mates of the covariance matrix verifies that the low level of correla- 
tion expected between each group of 9 or 4 simulations, introduced 
by splitting each fully independent simulation into sub-fields, does 
not impact significantly on our results. It also demonstrates that the 
Anderson] ((2003) inverse covariance de-biasing correction (equa- 
tion|16|, is sufficiently accurate for our analysis. The second con- 
servative analysis is to remove the angular scales where uncertainty 
in the accuracy of the non-linear correction to the power spectrum 
could bias our results. We select these angular scales by calculating 
a WMAP7 cosmology theoretical prediction for ^ (0) assuming 
two different non-linear corrections, where we boost and decrease 
the |Smithet^] l |2003) non-linear correction to the power spectrum 
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by ±7 per cent. Note that we chose the value of 7 per cent from the 
average error over the range of k scales tested in Eifler ( 201 1}. For 
angular scales where more than a 10 per cent difference is found in 
the expected signal, between these two different non-linear correc- 
tion regimes, we remove these scales from our analysis. As the 
statistic probes significantly smaller k scales compared to the 
statistic, at a fixed 8, we cut more £_ data than (see |Benjamin| 
|et al.|2013] for further discussion). For £+, our requirement for less 
than a 10 per cent deviation corresponds to the removal of data with 
< 3 arcmin for tomographic bin combinations including bins 1 
and 2. For , this corresponds to removing data with 6 < 30 ar- 
cmin for tomographic bin combinations including bins 1, 2, 3 and 
4, and data with 8 < 16 arcmin for tomographic bin combinations 
including bins 5 and 6. Applying these cuts in angular scale results 
in a final data vector of length p = 120. As the £± statistic is an 
integral over many k scales weighted by Jo and J4 Bessel func- 
tions, one cannot directly relate the limits we place on 0, to limits 
on k. We note, however, that as these cuts do preferentially remove 
the smallest physical k scales where the theoretical prediction to 
the power spectrum is expected to be most impacted by baryonic 
feedback effects. This conservative analysis to test the non-linear 
correction therefore also works as a mitigation strategy to avoid 
uncertain baryon feedback errors. For this conservative analysis we 
find no change in the best-fit measurement of <78(fl m /0.27) a , but a 
reduction in the constraining power by roughly 20 per cent (see the 
'Low 8 scales removed' row in table [2](. We also lose roughly 20 
per cent of the constraining power on the intrinsic alignment am- 
plitude A with this conservative analysis. As the best-fit value for 
(J8(r2m/0.27) a remains unchanged, we can conclude that the in- 
clusion of small-scale data does not introduce any significant bias 
in our results. Furthermore, as our focus for this analysis is the mit- 
igation of intrinsic galaxy alignments, which are most tightly con- 
strained by the low-redshift bins preferentially cut with this type 
of conservative analysis, the CFHTLenS results that follow include 
the full angular scale range shown in Figure [2] 

4.3 Joint Cosmological Parameter constraints 

We present joint cosmological parameter constraints from 
CFHTLenS combined with WMAP7, BOSS and Rl 1 for four cos- 
mological models testing flat and curved ACDM and uiCDM cos- 
mologies. Table[3]lists the best-fit 68 per cent confidence limits for 
our cosmological parameter set for the combination of CFHTLenS 
and WMAP7 (first line for each parameter), CFHTLenS, WMAP7 
and Rll (second line for each parameter) and for CFHTLenS, 
WMAP7, BOSS and Rll (third line for each parameter). For 
comparison the figures in this section also show constraints for 
WMAP7-only and WMAP7 combined with BOSS and Rll. We 
refer the reader to Komatsu et al.| <|20 1 1) and Anderson et al.| 
([20T2J for tabulated cosmological parameter constraints for the 
non-CFHTLenS combination of data sets shown, noting that we 
find good agreement with their tabulated constraints. We also re- 
fer the reader to Kilbinger et al. (2013) for CFHTLenS-only pa- 
rameter constraints for the curved and wCDM cosmological mod- 
els tested in this section. Whilst CFHTLenS currently represents 
the most cosmologically constraining weak lensing survey, it spans 
only 154 square degrees and is therefore not expected to have sig- 
nificant constraining power when considered alone. This is demon- 
strated in Figure [5] which compares parameter constraints in the 
erg — fi m plane for a flat ACDM cosmology. The wide constraints 
from CFHTLenS-only are shown in pink (note the inner 68 per 
cent confidence region was shown in pink in Figure|4j, in compar- 
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Figure 5. Flat ACDM joint parameter constraints (68 and 95 per cent 
confidence) on the amplitude of the matter power spectrum controlled by 
erg and the matter density parameter £7 m from CFHTLenS-only (pink), 
WMAP7-only (blue), BOSS combined with WMAP7 and Rl 1 (green), and 
CFHTLenS combined with BOSS, WMAP7 and Rl 1 (white). 

ison to the tight constraints from WMAP7-only (blue). The power 
of lensing, however arises from its ability to break degeneracies 
in this parameter space owing to the orthogonal degeneracy direc- 
tions. BOSS combined with WMAP7 and Rl 1 is shown green and 
when CFHTLenS is added in combination with BOSS, WMAP7 
and Rll (white) we find the combined confidence region decreases 
in area by nearly a factor of two. As we will show in this section, 
the tomographic lensing information presented in this analysis is 
therefore very powerful when used in combination with auxiliary 
data sets. 

The figures that follow in this section all compare constraints 
for different combinations of cosmological parameters and cosmo- 
logical models with the following colour-scheme: WMAP7-only 
(in blue), WMAP7 combined with CFHTLenS and Rll (in pink), 
WMAP7 combined with BOSS and Rll (in green) and all four 
data sets in combination (in white). Comparing the green contours 
with the pink contours allows the reader to gauge the relative power 
of BOSS and CFHTLenS when either survey is used in combina- 
tion with WMAP7 and Rll. Comparing the green contours with the 
white contours allows the reader to gauge the extra contribution that 
CFHTLenS makes to BOSS, Rl 1 and WMAP7 in breaking differ- 
ent parameter degeneracies and constraining cosmological param- 
eters. 

4.3.1 Constraints in the as — f2 m plane 

Figure|6]shows joint parameter constraints on the normalisation of 
the matter power spectrum erg and the matter density parameter 
Q m for four cosmological models: flat ACDM, flat toCDM, curved 
ACDM and curved mjCDM. The comparison of the results for the 
four cosmological models show the decreased WMAP7 sensitiv- 
ity to these two cosmological parameters when extra freedom in 
the cosmological model is introduced, such as dark energy wq, or 
curvature. We find slightly tighter constraints from CFHTLenS in 
combination with WMAP7 and Rll (pink), in comparison to BOSS 
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Table 3. Joint cosmological parameter constraints for four models, testing flat and curved ACDM and uiCDM cosmologies. The first line for each parameter 
lists the constraints from CFHTLenS combined with WMAP7 and Rl 1. The second line lists the constraints from CFHTLenS combined with WMAP7, BOSS 
and R 1 1 . Deduced parameters are indicated with * . 



Parameter flat ACDM flat toCDM curved ACDM curved ujCDM Data 
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in combination with WMAP7 and Rll (green). The 68 per cent 
confidence regions between these two survey combinations only 
marginally overlap, introducing a mild tension. The constraints are 
however consistent at the 95 per cent confidence level. For the mat- 
ter density parameter S7 m , the addition of BOSS data to the com- 
bined CFHTLenS, WMAP7, Rll analysis typically decreases the 
1<7 errors by ~ 20 per cent across all cosmologies. For the normal- 
isation of the matter power spectrum as, however, we find BOSS 
adds little to the constraining power of CFHTLenS with WMAP7 
and Rl 1 for the cosmological models tested. Furthermore, for a flat 



ACDM cosmology the constraint os = 0.799 ± 0.015 is almost 
entirely driven by CFHTLenS in combination with WMAP7 alone. 

4.3.2 Curved cosmological models 

We consider two curved cosmologies where the sum of the dif- 
ferent density components of the Universe is no longer limited 
to the critical density. Figure [7] shows joint parameter constraints 
on the curvature JIk and the matter density parameter f2 m for 
WMAP7-only (blue), BOSS combined with WMAP7 and Rll 
(green), CFHTLenS combined with WMAP7 and Rll (pink) and 
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Figure 6. Joint parameter constraints on the normalization of the matter power spectrum erg and the matter density parameter f2 m from WMAP7-only (blue), 
BOSS combined with WMAP7 and Rll (green), CFHTLenS combined with WMAP7 and Rll (pink) and CFHTLenS combined with BOSS, WMAP7 and 
Rll (white). 



CFHTLenS combined with BOSS, WMAP7 and Rll (white). In 
both the curved ACDM and curved iuCDM cosmology we find that 
the data are consistent with a flat Universe with £Ik ~ —0.004 ± 
0.004 (see Table[3]for exact numbers for the different cosmologies 
and data combinations). In this parameter space we find a factor 
of two improvement when Rll is included in combination with 
CFHTLenS and WMAP7. This is partly because when curvature is 
allowed the degeneracy direction of the CMB in the erg — S7 m plane 
changes such that the combination of lensing with the CMB be- 
comes less powerful. Little improvement is found in the constrain- 
ing power when BOSS is included in our parameter combination, 
but the mean Q,k changes by nearly 2a. 



where a 2a difference is found between the mean wq measured 
with lensing, WMAP7 and BOSS, with and without a prior on the 
hubble parameter h. For all cosmologies tested in this analysis, the 
constraints on h from CFHTLenS with WMAP7 are in good agree- 
ment with the Rl 1 measure of h = 0.738±0.024. Focussing on flat 
ACDM, in this analysis we find h = 0.717 ± 0.016 for CFHTLenS 
with WMAP7, in comparison to BOSS with WMAP7 who find a 
2a offset from Rl 1 with h = 0.684 ± 0.013. For wCDM cosmolo- 
gies, [Kilbinger et al.|(20T3] > and BOSS find even larger shifts away 
from the Rl 1 result, but at a lower significance, and it is this that 
causes the differences in the measurement of w , with and without 
the inclusion of a prior on h, that we do not find in this analysis. 



4.3.3 Constraints on dark energy 

Finally we turn to the constraints that can be placed on the dark en- 
ergy equation of state parameter wq in flat and curved cosmologies. 
Figure[8]shows joint parameter constraints in the w — fi m plane and 
also the w — Ok plane for a curved mCDM cosmology. As with the 
other parameter planes that we have commented upon in this Sec- 
tion, we again see the mild tension between BOSS and CFHTLenS 
and the power of including these surveys in addition to WMAP7 
data alone. For both the curved and flat wCDM cosmologies we 
find that w is consistent with a cosmological constant (see Table [3] 
for exact numbers for the different cosmologies and data combina- 
tions). As with our constraints on the curvature, we find very little 
improvement in the constraining power on w when BOSS is in- 
cluded in our parameter combination. We do, however, find excel- 
lent agreement with the combined probe constraints from |Ander-| 
|son et al.| ( pOT2l > when BOSS is combined with WMAP7, SDSS- 
LRG baryon acoustic oscillation constraints from Padmanabhan 
|et al.| p012) and Type la supernovae results from Conley et al. 
( |2011fr . With this combination of data sets, And erson et al.H2012) 
find w = —1.09 ± 0.08 for flat and curved wCDM models. Note 
that this was shown to be the only parameter where the addition of 
the supernova data to the BOSS and WMAP7 data impacted upon 
the analysis, decreasing the errors by a factor of ~ 2. This result is 
in agreement with our wCDM model constraints from CFHTLenS 
with WMAP7 and Rll, where we find w = -1.06 ± 0.08 (flat) 
and w = -1.04 ± 0.12 (curved). 

We find good agreement between the mean measurements 
when the different parameters sets are combined. This is in con- 
trast to the 2D weak lensing analysis of Kilbinger e t al.| J2013) 



5 THE INTRINSIC ALIGNMENT OF EARLY-TYPE AND 
LATE-TYPE GALAXIES 

As discussed in Section Q] there is clear evidence in the literature 
that the strength of the intrinsic alignment signal depends on galaxy 
type, with the most massive red galaxies exhibiting the strongest in- 
trinsic correlations ( Joachimi et al. 201 1 1. In this section we there- 
fore present separate tomographic analyses of an early-type and 
late-type galaxy sample, selected using the measured best-fit pho- 
tometric type Tbpz- This classification type ranges from 1 to 6 
and represents the best-fitting spectral energy distribution to each 
galaxy's photometry (see Hildebr andt et al.|2012| for more details). 
We follow Simon et al.]l |2013| > by selecting late-type spiral galaxies 
with Tbpz > 2.0, roughly 80 per cent of the galaxy catalogue used 
in the main analysis. The remaining 20 per cent are classified as 
early-type galaxies. Each sample is split into 6 tomographic bins, 
using the redshift selection given in Table[T] and the redshift distri- 
bution determined from the sum of the P(z). Covariance matrices 
were determined for each galaxy sample using the method outlined 
in Section [373] but mapping only the relevant galaxy sample on to 
the N-body lensing simulations. We found no evidence for a signifi- 
cant difference in a e for the two samples. Figure[9]shows the result- 
ing compressed tomographic measurements made with early-type 
galaxies (circles) and late-type galaxies (crosses). The data com- 
pression uses the visualization method described in Section [4~l] 
modified slightly such that the free amplitude parameter a v is fit- 
ted simultaneously to both the and £_ measured from the data. 
This simultaneous fit is justified as the a 13 values measured for 
£+ and £_ independently are fully consistent. The resulting best- 
fitting amplitude a lJ is shown, multiplied by the fiducial model at 
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Figure 8. Joint parameter constraints on the dark energy equation of state parameter wo and the matter density parameter f2 m , and curvature parameter f^K 
for a curved wCDM cosmology from WMAP7-only (blue), BOSS combined with WMAP7 and Rl 1 (green), CFHTLenS combined with WMAP7 and Rl 1 
(pink) and CFHTLenS combined with BOSS, WMAP7 and Rl 1 (white). 



6 = 1 arcmin for . With only 20 per cent of the data contained 
in the early-type sample, it is unsurprising that the measured signal 
to noise is significantly weaker than for the late-type sample which 
are well fit by the fiducial GG-only model, shown dotted. We can, 
however, optimise the measurement of the intrinsic alignment sig- 
nal from early-type galaxies, to get a clearer picture, if we assume 
the II contribution to cross-correlated bins is small in comparison 
to the GI signal. If this is the case, we can decrease the noise on 
the GI measurement by using the full galaxy sample as background 
galaxies to correlate with the early-type galaxies in the foreground 
bin. The result of this optimised analysis is shown, in compressed 
tomographic data form, in Figure [TU] The open circles show the 



tomographic signal measured in the auto-correlated redshift bins 
between early-type galaxies (these auto-correlation bins are also 
shown in Figure|9|. The closed symbols show the tomographic sig- 
nal in the cross-correlated redshift bins where early-type galaxies 
populate the foreground bin and the full galaxy sample populates 
the background higher redshift bin. The data can be compared to 
the fiducial GG-only model, shown dotted. What is interesting to 
note from this Figure is that at low redshifts, where the intrinsic 
alignment signal is expected to be the most prominent, the auto- 
correlated bins tend to lie above the GG-only model. We expect 
this from the II term. For the cross-correlated bins, however, the 
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Figure 9. Compressed CFHTLenS tomographic data for two galaxy sam- 
ples; early-type (circles) and late-type (cross) galaxies. As in Figure]!] each 
point represents a different tomographic bin combination ij as indicated by 
Zpeak' me peak redshift of the lensing efficiency for that bin. The measured 
best-fitting amplitude a IJ of the data for each galaxy type, multiplied by the 
fiducial model at 8 = 1 arcmin for is shown. The error bars show the 
lcr constraints on the fit. The data can be compared to the fiducial GG-only 
model, shown dotted. 

measured signal tends to lie below the GG-only model. We expect 
this from the GI term. 

FigureJTTJcombines the CFHTLenS data split by galaxy type, 
and our optimised early-type galaxy tomography analysis, with 
auxiliary data from WMAP7, BOSS and Rl 1 to constrain the am- 
plitude of the intrinsic alignment model A. Assuming a flat ACDM 
model, the resulting 68 per cent and 95 per cent confidence limits 
on A and the matter density parameter f2 m can be compared^ In 
the left panel we show constraints from the two galaxy samples 
split by SED type. The early-type galaxy constraints are shown 
in red and the late-type galaxy constraints are shown in blue. In 
the right panel, constraints are shown for the full galaxy sample 
in purple and the optimised early-type intrinsic alignment analysis 
in pink. The marginalised 68 per cent confidence errors on A, from 
the combination of CFHTLenS data with WMAP7, BOSS and Rl 1, 
for the four different measurements are 

Alste = 0.18±S:1I , (17) 
Nearly — 5.15^ 2 .32 1 (18) 

A^ ly = 4.26±\f 9 , (19) 

Asa = -0.48t ( o;s7 ■ (20) 

We find the intrinsic alignment amplitude of the late-type sample 
is consistent with zero. In contrast, the amplitude of the intrinsic 

5 Note that the constraints on cosmological parameters other than A are 
consistent between the early-type and late-type analysis, and that both sets 
of parameter constraints, with the exception of A, are consistent with the 
full galaxy sample analysis reported in table[51 



Figure 10. Compressed CFHTLenS tomographic data for an optimised 
early-type galaxy intrinsic alignment measurement with auto-correlated 
redshift bins containing only early-type galaxies (circles) and cross- 
correlation redshift bins containing early-type galaxies in the low redshift 
bin and all galaxy types in the high redshift bin (filled). Different tomo- 
graphic bin combinations ij are indicated by z pca k> the peak redshift of the 
lensing efficiency for that bin. The best-fitting amplitude a %3 of the data 
relative to a fixed fiducial GG-only cosmology model is shown, multiplied 
by the fiducial model at 8 = 1 arcmin for £+. The error bars show the la 
constraints on the fit. The data can be compared to the fiducial GG-only 
model, shown dotted. 

alignment model for the early-type sample is detected to be non- 
zero with close to 2a confidence. When we consider the optimised 
analysis, we find an even stronger detection, with an intrinsic align- 
ment amplitude of A = for early-type galaxies ruled out with 
3cr confidence. The optimised early-type analysis should be con- 
sidered with some caution, however, as the tomographic redshift 
bins do overlap and as such a small fraction of late-type with early- 
type II correlation will be included in the measurement. The mea- 
surement of A C ariy should therefore be considered as our cleanest 
measurement of the early-type galaxy intrinsic amplitude with the 
optimised A°^ tl analysis providing us with the strongest evidence 
for intrinsic galaxy alignments between early-type galaxies. 

5.1 Discussion 

Our constraints show the same broad findings as other studies; in- 
trinsic alignments are dependent on galaxy type. As previous stud- 
ies have focused on specific galaxy samples at fixed redshifts, how- 
ever, it is difficult to compare our constraints directly. With that 
caveat we can, however, comment on literature results from galaxy 
samples that are the most comparable. Our late-type sample is most 
similar in its properties to the blue galaxies from the WiggleZ sur- 
vey analysed in Mandelbaum et al. ( 20 1 1 1. Their null detection is in 
agreement with our late-type galaxy results. Our early-type sample 
is most similar in terms of luminosity and redshift to the MegaZ- 
LRG sample analysed in Joachimi et al. (201 1 1. The best-fit values 
4 < A < 6 for a range of different types of LRG galaxy selection 
with an error of ~ 1, are in very good agreement with our early- 
type galaxy results. 

For the full galaxy sample, there is an indication that negative 



© 0000 RAS, MNRAS 000, 000-000 



18 C. Heymans et al. 







| Optimised Early-type ■ 
[| Full sample 

* '*' "\ >. 
/ \ \ 

' i \ 
i j i 

\ i f } 

\ \ if 

\ N / / 


V^^^'^ [J Early- type 
□ Late-type 







0.25 



0.3 



0.35 



0.25 



0.3 



0.35 



Figure 11. Joint parameter constraints on the amplitude of the intrinsic alignment model A and the matter density parameter Q m from CFHTLenS combined 
with WMAP7, BOSS and Rl 1. In the left panel the constraints can be compared between two galaxy samples split by SED type, (early-type in red and late-type 
in blue). In the right panel we present constraints from a optimised analysis to enhance the measurement of the intrinsic alignment amplitude of early-type 
galaxies (pink). The full sample, combining early and late-type galaxies, produces an intrinsic alignment signal that is consistent with zero (shown purple). A 
flat ACDM cosmology is assumed. 



values of A are preferred. For flat cosmologies, A is negative at 
the 1.4a level when the CFHTLenS data are combined only with 
WMAP7 and Rll (see table [5] for constraints on A for the full 
galaxy sample for different cosmologies and data combinations). 
Whilst we emphasize that this result is not statistically significant 
it is however worth commenting on what this finding could mean. 
In the conventional intrinsic alignment model the GI signal is neg- 
ative and scales with A. The II signal is positive and scales with 
A 2 . Finding A < 0, however, implies the data prefer a GI+II sig- 
nal that is more positive than the conventional model would pre- 
dict. This suggests that future surveys with lower statistical errors, 
should aim to fit independent amplitudes to the GI and II signals as 
the interplay between the two effects may be more complex than 
the linear tidal field alignment model suggests. 

It is also interesting to comment on the decrease in the ampli- 
tude of the best-fit intrinsic alignment signal when early and late 
type galaxies are combined. If detected in future surveys at higher 
significance, this would indicate a complex interplay between the 
two galaxy types. It has long been thought that the reason for the 
difference between the intrinsic alignments of early and late type 
galaxies lies in the different mechanisms at play during galaxy for- 
mation. The intrinsic alignment model we use in this analysis is 
based on linear theory. A more traditional galaxy formation sce- 
nario for late-type galaxies, however, is tidal-torque theory where 
it is the angular momentum of the dark matter field that induces 
galaxy spin and hence intrinsic galaxy alignments (see [Schafer 



2009 and references therein). The simple hypothesis, presented in 



Heym ans et al.| l (2"006[ >, is that the intrinsic alignment of early-type 



galaxies is a result of ellipticity deformations due to the linear tidal 
field, in contrast to late-type galaxies whose alignment results from 
angular momentum induced ellipticity alignments ( van den Bosch 
|et al.|2 002 ). This hypothesis is in good agreement with recent ob- 
servations of galaxy-type dependence in the intrinsic alignment sig- 
nal, as halo angular momentum is proportional to the square of the 



tidal shear, and the induced galaxy alignments therefore correlate 
over much shorter ranges compared to alignments directly caused 
by the linear tidal shear (Ca telan et al.|200T| >. 

In addition to the linear model used throughout this paper, |HTj 
rata & Seljak (2004) also investigate the GI signal expected from 
an intrinsic alignment model where the galaxy ellipticity is pro- 
portional to the square of the tidal field. In this case the GI sig- 
nal is expected to be zero. As our galaxy sample is dominated by 
late-type galaxies, the majority of correlated galaxy pairs in our 
analysis from different redshift bins will include a late-type fore- 
ground galaxy. Combining the findings of Hirata & Seljak {2004| > 
with our simple hypothesis that late-type intrinsic galaxy align- 
ment is caused by halo angular momentum induced alignments, 
leads to an expected zero GI measurement on average. In auto- 
correlated tomographic bins however, the stronger galaxy cluster- 
ing of early-type galaxies will mean that at small angular scales, 
there is a higher proportion of close early-type galaxy pairs in the 
measurement, compared to the numbers of early-type and late-type 
foreground galaxies that contribute to the GI signal. This therefore 
boosts the true II signal in auto-correlated bins over the amplitude 
that would be predicted from Gl-only constraints from a mixed 
galaxy population. 

The linear tidal field alignment model used in this analysis 
could compensate for these different galaxy-type contributions to 
the II and GI signal by favoring a small but negative value for A. 
In this case the GI signal in the cross-correlation bins is positive 
but sufficiently weak to provide a reasonable fit to the GI=0 model 
signal expected from the dominant late-type galaxy population. In 
the auto-correlated bins, the additional true positive II signal from 
the clustered early-type galaxies is then represented in the model 
fit, not by the model II signal, but the positive GI signal. If A was 
positive and less than unity, there would still be a reasonable weak 
but now negative fit to the GI=0 model in the cross-correlation bins. 
In the auto-correlated bins, however, there would not be sufficient 
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signal in the combined II+GI model to represent the extra II power 
arising from the clustered early-type sample. 

Based on this discussion, we can conclude that our constraints 
for the full sample favoring a slightly negative value for A fits our 
simple hypothesis that early-type galaxy alignment results from the 
linear tidal field and late-type galaxy alignment results from angu- 
lar momentum induced correlations. The next generation of weak 
lensing surveys will have the statistical power to test this hypothesis 
further. 



6 CONCLUSIONS 

The Canada-France-Hawaii Telescope Lensing Survey, 
CFHTLenS, represents the current state-of-the-art in cosmo- 
logical weak lensing data analysis, from the applied weak lensing 
optimised data reduction, shear and photometric redshift measure- 
ment methods, through to the robust systematic error analysis and 
error quantification of the resulting shear and redshift catalogue. 
Spanning 154 square degrees, CFHTLenS is currently the largest 
deep weak lensing survey in existence permitting the tightest 
cosmological constraints from weak gravitational lensing. In this 
article we present the first multi-redshift bin, or tomographic, weak 
lensing analysis to mitigate the contamination to the measured 
two-point shear correlation function through the simultaneous fit 
of a cosmological model with an intrinsic galaxy alignment model. 
Combining the tomographic CFHTLenS data with auxiliary 
cosmological probes; the cosmic microwave background with 
data from WMAP7, baryon acoustic oscillations with data from 
BOSS, and a prior on the Hubble constant from the HST distance 
ladder, we have improved constraints on a range of cosmological 
parameters for a standard flat ACDM model, in addition to curved 
and dark energy models. We constrain the amplitude of the matter 
power spectrum trg = 0.799 ± 0.015 and the matter density 
parameter tt m = 0.271 ± 0.010 for a flat ACDM cosmology. For 
a flat iuCDM cosmology we constrain the dark energy equation 
of state parameter w — —1.02 ± 0.09. In general we find tighter 
constraints from the combination of CFHTLenS with WMAP7 
and Rl 1 than from BOSS with WMAP7 and Rl 1 and we find that 
the addition of BOSS to CFHTLenS with WMAP7 and Rll only 
significantly improves constraints on the matter density parameter 
S7 m , for all cosmologies tested. Constraints on the other parameters 
are only shown to significantly improve when a curved u>CDM 
model is considered. Finding consistent results, however, between 
these two very different probes of cosmology suggests a bright 
future for studies of the 'Dark Universe' with weak lensing and 
baryon acoustic oscillations. 

Tomographic weak lensing has long been recognized as a 
powerful tool to constrain dark energy by detecting the influence 
dark energy has on the growth of structure in addition to the 
distance-redshift relationship. One astrophysical source of uncer- 
tainty that mimics cosmological weak lensing is the intrinsic align- 
ment of neighbouring galaxies. This phenomenon unfortunately re- 
duces the overall constraining power of tomographic weak lensing 
analyses as, to ensure the cosmological constraints are unbiased, 
the contamination from intrinsic alignments must be considered. In 
this analysis we have assumed a simple one-parameter model that 
scales the amplitude of the II and GI intrinsic alignment contami- 
nation expected from a linear tidal field alignment model of galaxy 
shape correlation jCatelan et al.|2001||Hirata & Seljak|2004||Bri-| 
|dle & King 20071. Our results are consistent with there being zero 
intrinsic alignment between galaxies in our late-type sample, with 



A ~ 0.2 ± 0.8. The best-fit amplitude is, however, within la of the 
fiducial model amplitude A = 1. This fiducial model is commonly 
assumed for parameter forecasts, and is based on the amplitude 
of low-redshift galaxy intrinsic ellipticity correlations measured by 
|Brown et al.| ( [2002"| >. For the 20 per cent of galaxies in our sample 
whose 5-band photometry is best-fit by an early-type SED, how- 
ever, we detect a non-zero intrinsic alignment signal, A ~ 5 ± 2, 
roughly five times the fiducial model amplitude. This is in agree- 
ment with previous observations of a galaxy-type dependence, us- 
ing a very different methodology ( Hirata et al. 2007 ; Joachimi et al. 
|2011[ [Mandelba um et al.||2011) . Our results therefore add to the 
increasing body of independent observations that point towards a 
scenario where the galaxy formation and evolution mechanisms, 
which determine galaxy shape, differ for early and late-type galax- 
ies. For the combined galaxy sample, we find that the net effect of 
the two galaxy types produces an intrinsic alignment signal that is 
consistent with zero, with A ~ —0.5 ± 0.8, (see Section[4]for the 
exact constraints for the different cosmological models). 

One difficulty, for the analysis of finely-binned tomographic 
correlation functions which are optimised for the simultaneous 
analysis of cosmological models and intrinsic alignment models, 
is the estimation of an accurate and invertible covariance matrix for 
the large data vector. In this paper we discuss the different options 
for covariance matrix estimation and the noise biases that can arise 
in the matrix inversion. We chose to use N-body lensing simula- 
tions as the basis for our covariance matrix estimation in order to 
correctly account for non-Gaussianity on small angular scales. This 
is in contrast to many earlier studies which often resorted to a fitting 
function correction. This choice, however, limits the largest angular 
scales that we can analyse. The finite simulated box size truncates 
the large-scale modes, reducing the large scale power probed by 
each simulated line of sight. In addition, as high-resolution N-body 
lensing simulations are expensive to create, we are also limited by 
the total number of simulated independent lines of sight available. 
This finite number limits the number of data points for which we 
can invert the simulation estimated covariance matrix without bias- 
ing our results, or erroneously increasing the area of the resulting 
confidence regions. A priority of future surveys must therefore be to 
ensure the availability of a large volume of N-body simulations for 
covariance matrix estimation. Alternative methods should also be 
developed such as the simulation-analytic covariance matrix graft- 
ing method described in Kilbinger et al. ( 20 1 3 1, in addition to meth- 
ods which optimally combine simulation and data-based bootstrap 
or jackknife estimates. 

Three major new weak lensing survey^] are currently in 
progress, or will commence soon, and as such this is an exciting 
time for the study of the 'Dark Universe'. These surveys will im- 
age more than ten times the area surveyed by CFHTLenS, chart- 
ing a sufficient volume to address the issue of astrophysical bias 
arising from intrinsic galaxy alignments in greater detail. By self- 
calibrating the intrinsic alignment model with the additional in- 
formation from galaxy clustering (Bernstein 2009; Zhang 2010; 
Joachimi & Bridle 2010), the marginalisation over many intrinsic 
alignment nuisance parameters will become feasible. Another vi- 
able mitigation strategy is to develop new statistics which are less 
sensitive to intrinsic alignment contamination. In addition, com- 
plementary spectroscopic or highly accurate photometric redshift 
observations could be exploited to place observational constraints 



6 KiDS: kids.strw.leidenuniv.nl, DES: www.darkenergysurvey.org and 
HSC: www.subarutelescope.org/Projects/HSC 
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on different analytical and hydro-dynamical simulation models of 
galaxy shape correlations. This effort, in parallel to large-area sur- 
vey acquisition, will then allow us retain as much power as possible 
from weak gravitational lensing as a cosmological probe. 
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